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WIND-TUNNEL LIFT INTERFERENCE ON SWEPTBACK WINGS IN 
RECTANGULAR TEST SECTIONS WITH SLOTTED SIDE WALLS 

By Ray H. Wright and James D. Keller 
Langley Research Center 

SUMMARY 

A theory is presented for the boundary-induced upwash interference on a sweptback 
wing mounted in a rectangular wind-tunnel test section with closed top and bottom walls 
and slotted side walls. The interference factor can be computed at any point in the test 
section. In an example, the interference with slotted side walls is compared with that 
with slotted top and bottom walls. 


INTRODUCTION 

In reference 1 a theory was derived for the subsonic wind-tunnel-boundary lift inter- 
ference on sweptback wings mounted in the horizontal center plane of rectangular test sec- 
tions with slots in the top and bottom walls, where the slots were used to reduce the inter- 
ference effects of the boundaries. However, it is sometimes convenient in making yaw 
tests to rotate the model through a right angle about the longitudinal axis of the test sec- 
tion. If the direction of positive lift is designated "up," the slots are then in the sides of 
the test section rather than in the top and bottom walls. In general, the change in con- 
figuration may be expected to produce a change in both magnitude and distribution of the 
boundary-induced upwash. The theory of reference 1 has been modified to apply to a 
sweptback wing in a test section with slotted side walls and solid top and bottom walls. 

The theory herein presented has been left in more general form than that of reference 1 
in that the upwash interference factor can be computed for any point in the test section 
and the model may be mounted anywhere in the test section, as long as the lift may be 
taken as (positive or negative) up. 


SYMBOLS 

a segment of airfoil span included in a discrete -point representation 


B(w,f) ,C(u;,f) 


parameters in solution of transformed Laplace equation 



b width of test section 

C cross-sectional area of test section 

Cl lift coefficient of model 

d distance between two adjacent slots 

F exponential Fourier transform of cp on x - Xj and 

F^' exponential Fourier transform of on z - z^ 

Fj' exponential Fourier transform on x - x^ and z - z^ 

potential corresponding to F^’ 

f variable of transformation on z - z^ 

h height of test section 

i = 

j summation index 

Kj modified Bessel function of the second kind 

k odd summation index 

L total lift of model 

AL element of lift 

l restriction constant of slotted walls 

M Mach number 

n . . . N even summation indices 


p = ho> 

Qij 


function of q defined by equation (26) 


Z - Zj 


of an interference 
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q = hf 


S 

V 


X,Y,Z 

x,y,z 

xijYi.Zj 

a 

T 

6 

6(o>) 

(p 

n 

<Pd 

<Pd 

*A 

O) 


ratio of slot width to distance between slot centers 

area on which is based 

tunnel stream velocity 

upwash velocity 

Cartesian coordinate axes 

Cartesian coordinates 

coordinates of image of a lifting element 

a positive real parameter 

circulation 

upwash interference factor 
Dirac delta 

coordinates of a lifting element 
velocity potential 
interference potential 

velocity potential of a semi-infinite line doublet 
modified velocity potential defined by equation (3) 
exponential Fourier transform of cp^ on x - x^ 
variable of transformation on x - Xj 
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Subscripts: 


i,j position identifications 


ANALYSIS 


The lifting wing is represented by a distribution of semi-infinite doublet lines 
starting at discrete points of lift application on the sweptback wing and extending down- 
stream toward infinity. The total interference is the sum of the interferences due to the 
interaction of the boundaries on the individual doublet lines. The doublets are oriented 
as indicated in figure 1, which also shows the coordinate system relative to the test sec- 
tion of height h and width b. With the origin at the center of the test section, the 
coordinate x is in the stream direction, y is normal to the slotted side walls, and 
z is vertical, positive upward. Positive lift is taken in the positive z-direction. The 
boundary-induced upwash velocity, positive in the positive z-direction, due to a single 
lifting element located at some point in the test section is now derived. 



Figure 1.- Schematic drawing of test section and doublet 
configuration showing coordinate designation. 


In a development similar to that in reference 1, the velocity potential at (x,y,z) 
due to a lifting semi-infinite line doublet starting at (x^y^z-jJ is found to be 


<?d 


JTa 
47 r 


z - z^ 

(y - yi) 2 + (z 



X - Xj 

]J(x - Xi) 2 + (y - yi) 2 + (z - Z X ) 2 


( 1 ) 


4 



where r is circulation and a is the segment of airfoil span over which the element 
of lift is carried. 

In order to facilitate the required integral transforms, let a modified potential be 


„ , _ ra 
^d "27 


z - z 1 


47r (y - yj) 2 + (z - zj) 2 


1 de ' a \l( x ~ x i) 2+ ( y - y if + ( z ' z i y 

1_ “ 8(X - x x ) 


( 2 ) 


where a > 0. Then 


U = %% < < 3 > 

and as in reference 1, the exponential integral transform of cp ^ on x - x^ with vari- 
able of transformation u> is 


n d'^>y» z ) = 


Z - Z-, 


(y - y x ) 2 +(z - 


[ 7r6(co) 


ia ) |/(y - yQ 2 + ( z - z i) ; 


Vco 2 


+ a' 


X K X 



(4) 


for a > 0 and y^y - y^) 2 + (z - z^ 2 > 0 where Kj is the modified Bessel function 
of the second kind. 


From reference 2 (formula 15, p. 65 and formula 43, p. 112), the exponential trans- 
form of O d ' on z - Zj with variable of transformation f is 


F d T (co,y,f) = 


i iifi - y-yi V w2 + f2 + a2 

. x f - y-yi f wf e 1 11 

ijr6(oj) -t— r e 1 ' + -* - = — 

l f l u 2 + a 2 y w 2 + f 2 + a 2 


(5) 


With transformations on x - x^ and z - z-^, the Laplace equation 
= 0 transforms to 

9x2 9y2 9 Z 2 

= (a> 2 + f 2 )F 

9y 2 

A solution of this equation is 

Fj'(w,y,f) = B(co,f) sinh(y)/a> 2 + f 2 )+ C(w,f) cosh(y\/a> 2 + f 2 ) 


( 6 ) 

(7) 
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As in reference 1, but with y taking the place of z since the side walls rather than 
the top and bottom walls are slotted, the equivalent homogeneous wall boundary condition 
(originally derived in ref. 3) is given by 


where the positive and negative signs apply at y = - and y = - -, respectively. The 
restriction constant l is given by 


( 8 ) 


d 


Trr, 


l = - 10 g e CSC — 


(9) 


where d is the distance between the centers of two adjacent slots and r Q is the ratio 
of slot width to the distance d. The boundary condition (8) transforms to 


F ±1 0 

ay 


( 10 ) 


With Fjj’ the transform of the modified potential of the doublet line, the trans- 
form Fj' of the corresponding interference potential is to be so determined that 
F = F(j' + Fj' satisfies the boundary condition (10). For this purpose the partial deriva- 
tives of F<j' and F^' at the slotted boundaries are required. 

From equations (5) and (7) the partial derivatives with respect to y for positive 
y - y^ are found to be 


f 8 *d\ _ ra 

i ,6( t c)fe'( y ' yi ) lf ' + e-( y - y >>' Uf2+ “ 2 

\ 9 y /y-y^o 2 

_ + ofl 


( 11 ) 


and 


8F,’ 


9y 


— = B(w,f) \/c o 2 + f 2 cosh y^a) 2 + f 2 + C(co,f) ^u> 2 + f 2 sinh y^L 2 + f 2 (12) 


Similarly, for negative y - yj 


fd 

^ 9y 


/y- yi <0 


Ta 

2 


i7r5(ct>) fe 


(y-yj)|f| 


cot 




(y-yi))/ a,2 +f 2 +« 2 


(13) 


and equation (12) again applies. Since y. lies between the boundaries, equation (11) 

i- 

applies at the boundary y = - and equation (13) applies at the boundary y = - H. Inser- 

tion of F^' + Fj for F in equation (10) at the two boundaries y = £ and y - _ £ 
gives, respectively, 
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Addition and subtraction of equations (14) and (15) yield, respectively, 



Solution of equations (16) and (17) for B(w,f) and C(w,f), respectively, and sub- 
stitution of the expressions so obtained into equation (7) gives the transform of the modi- 
fied interference potential as 
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F-, f = — 


i7T6(co) e ^ \ sinh^y 1 |fQ^-j£-j - lij + — — - - sinh^ \lco 2 + f 2 + ^ - lj sinh(y\|od 2 .+ f 2 ) 


-§\jw 2 +f 2 +oi 2 


sinh^w 2 + f 2 ^ + l fa 2 + f 2 cosh^^w 2 + f 2 j 


-k|f| / , -^CJ 2 +f 2 +Q! 2 


i7r6(w) e 2 cosh(y 1 |f|)^ - Ifj ■ 


wfe 2 


2 , cosh( yi \M + f 2 + a 2 Y— L - z) 

“ 2 + “ 2 V / \t 2 + l 2 + a 2 21 


cosh (y Vo) 2 + f 2 ) 


( 18 ) 


cosh^||u>2 + f2 j+ i^oj^ + f 2 sinh^^aj2 + f 2 ^ 


The interference potential due to the slotted side walls in the presence of the poten- 
tial <?£ of the lifting doublet line is now obtained by taking the limit as a approaches 
zero of the inverse transform of Fj_'. 


<?i = Um r <*<*> r 

o?— 0 8tt z 


oo H 

e 

-oo 


i[(X-Xi) W ( Z - Zl )f] a 


(19) 


and the upwash velocity is therefore 


dcp 1 

dz 


r _ r°° C 00 iMX-xi)o>+fz-zi)f| 

l im Ja dw e Lt V ( ^JifF/ 

a -0 8 71^ ^-°o ^-oo 


( 20 ) 


By taking the limit in equation (20) , integrating the term in F^’ containing S(o>), and 
noting the evenness in f, the upwash velocity v is found to be 


v = 


dcp l 

dz 


ra 

87T 2 


Ulfl- A"2 |f| 


[" sinh^yjl f | j sinh(y|f|) 

cosh^y ]f ) cosh(yf) ~j 

sinh 

jl f l) + z| f l cosh 


cosh 

(|f) + i|f|sinh(||f|) 


Y 1 

\Va>2 + f2 / 


sinjjx - x^cJj 

sinh^y^w 2 + f 2 jsinh(yta 2 + f 2 ) 

a; 

sinh^a> 2 + f 2 J + Z P + f 2 cosh^\/a> 2 + f 2 j 


cosh 

( y i^ 

. 2 + f z ) 

cosh 

j)P 

! + 

f 2 ) 


cosh 


+ f 2 

^ + l ^2 + f 2 8^ 


k 2 " 

H 2 ) 


dco > cos 


[(z - z,)f] 


df 


( 21 ) 
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The upwash interference factor is defined as 


5 - < 22 > 

where 

v upwash velocity 

V tunnel stream velocity 

C cross-sectional area of test section 

Cl lift coefficient of model 

S area on which Cl is based 

As in reference 1, 


C Ta _ C AL / 2 o\ 

SVC l8 7t2 16tt2 L 

where AL is the element of lift represented by Ta and L is the total lift of the 
model. Then with C = hb and with change of variables to p = hw and q = hf, the con- 
tribution to the upwash interference factor at a point (x,y,z) due to the slotted-wall influ- 
ence associated with a lifting element at (xj,yj,z^ is found from equations (21), (22), 
and (23) to be 



-> OO 

o L 

sinh (ir q ) sinh (h q ) 

cosh (~ qj cosh^ 

8»HL 

_ sinh (l S q ) - E q cosh (i S ' 

q ) cosh (l K q ) + h q sinh (l h q ) 


^oo 

(l q2 q 2 ^ 

-~'|p 2 -hi 2 sin 


sinh^^p 2 + q 2 j sinh^j^p 2 + q 2 ^ 

J„\ h 

! p 

sinh0 j^p 2 + q 2 ) + ^p 2 + q 2 cosh^l hfa 2 + q 2 j 



j 
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I 


In equation (24) lim 
P“-0 


fe-£>] _ £ .fi 


The limits as q or q and p approach 


zero of all other expressions with zeros in the denominator are zero. The integrals 
should be convergent in the neighborhood of p = 10 and q = 20. 


Let the original lifting element be located at (£, 77 ,^); then the element and its images 
in the top and bottom walls are alternately positive and negative and are located at 
(£,? 7 ,kh - £) and at (|,rj,nh + ?), where 


k = ±1, ±3, ±5 . . . 

and 

n = 0, ±2, ±4, ±6 . . . 


To every such element there corresponds a contribution A5 to the upwash interference 
factor arising from the slotted side walls and having the form of equation (24). Further 
contributions to the downwash interference factor are made by the images in the solid 
top and bottom walls. From differentiation of equation (1) the upwash velocity produced 
at (x,y,z) by a lifting element located at is 



By equations (22) and (23), the corresponding contribution to the upwash interference fac- 
tor is 
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Let (A5)jj be the contribution to the upwash interference factor at a 
point corresponding to a lifting element located at a point Then by 

summing all the contributions of the form of equation (24) due to the slotted side wall? 
and all contributions of the form of equation (25) due to the images in the top and bottom 
walls, the total contribution to the upwash interference factor at (xj,yj,z^ due to the lifting 
element at is obtained. For convenience in the summation, let 



Then 





where 



k = ±1, ±3, ±5 . . . 
n = 0, dt2, ±4, ±6 . . . 


11 



The total upwash interference factor at due to all lifting elements in the pres- 

ence of the tunnel boundaries is 

64 = ^ 6 )^ (28) 

i 

If the lift distribution of the model lies symmetric with respect to the XZ-plane, 

the terms of containing sinh^Aqj and sinh^j^p^ + can be ignored, since 

their effects will add to zero in the summation of equation (28). Also with this type of 
symmetry, values of y^ can be restricted to positive or to negative values, since the 
interference is likewise symmetrical relative to the XZ-plane. 

If the lift distribution may be considered to lie in the XY-plane, the summations on 
n and k can be combined, because is zero. 

SAMPLE CALCULATION 


Let a wing swept back 35° and spanning 0.7 of the test section width be mounted at 
the center of a square, slotted wind-tunnel test section with four symmetrically spaced 
slots in each of the side walls and with solid top and bottom walls as shown in figure 2 . 


z 



Figure 2.- Schematic drawing of cross section 
of wing and test section. For sample calcu- 
lation, sum of slot widths in boundary of 
height h is O.C>6h. 


Let the slots occupy 6 percent of the area of each slotted boundary. Then by application 
of equation (9) , the restriction constant l is given by 


l 

h 


- is Iog e csc nr = w log < 


CSC 


0.06ir 

2 


_ 1 b i 

4 it h 1 g e 


csc(0.037r) ~ 0.18 


Let the lifting wing be represented by lifting elements located at points Pj, P 2 , . . P 10 

on lines of 35° sweep as shown in figure 3. 
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Side walls of test section 


i 



Figure 3-- Representation of sweptback wing 
for sample calculation. 


The coordinates of these points and the lift distribution assumed are given in the following 
table: 


Point 

Cj x i 

h ° r h 

^ nr * 
h h 

Sj 

h 

z i 

° r h 

m, 

P 1 

0.0246 

0.0351 


0 

0.1342 

P2 

.0738 

.1054 


0 

.1334 

P3 

.1229 

.1756 


0 

.1118 


.1721 

.2458 


0 

.0769 

P5 

.2212 

.3160 


0 

.0437 

P6 

.0246 

-.0351 


0 

.1342 

P? 

.0738 

-.1054 


0 

.1334 

P 8 

.1229 

-.1756 


0 

.1118 

P0 

.1721 

-.2458 


0 

.0769 

P 10 

.2212 

-.3160 

. 

0 

.0437 


\ The values in this table are the same as those used for the sample calculation in refer- 

ence 1. The lift distribution given here approximates that on a model which was of inter- 
est when the original calculation was performed. 
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With substitution of £ = 1 and £ = 0.18 

h h 


and with use of the coordinates 


Yi Zi /at\ 

-J-, — , and — and of the lift distribution (-j— j . 


h’ h’ h’ h 
interference factor 


h’ h’ 

given in the table, the upwash 

(A 6)^ at any point Pj corresponding to a lifting element at any 
point Pj can be computed by use of equation (27) where p and q must be carried to 
values large enough to insure convergence of the integrals with infinite upper limits and 
k and n must be carried to values large enough to insure convergence of the summa- 
tions. The total upwash interference factor at point is the sum of contributions from 
all points Pj and is given by equation (28) as 


5i = y (i 5 ) i3 


Lj 

j 


Since the wing is symmetrical about its midspan, it was necessary to compute for 
points on only one side of the midspan. The points Pj, P2, P3, P4, and P5 were chosen. 

Computation of was made also for the points Pq = ^ = (0,0) and 

Pjj = ^ = (0.2451, 0.350). (Note that Po and Pjj are i-points, but not j-points, 

whereas all other points are both i-points and j-points.) The calculation was performed 
on a Control Data 6600 digital computer. 

The calculated values of 6i are shown in figure 4 as a function of spanwise posi- 
9 Yi 

tion jpj on units of the semispan. For comparison, the upwash interference factor 
for the same wing and lift distribution in the same test section, but rotated 90° about the 
streamwise X-axis (or, otherwise, slots in top and bottom walls), is also shown from 
reference 1. 



Figure 4.- Upwash interference factor for wing 
swept back 35° and- spanning 0.7 of width of 
a square test section with four slots in 
each of its slotted walls, slot opening 
being 6 percent of each slotted wall. 
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DISCUSSION 


The effects of rotating the lift distribution given in the sample calculation through 
a right angle about the X-axis are seen to be: 

1. With the slots on the sides instead of in the top and bottom walls, the sign of 

the lift interference was reversed over most of the span. The average interference over 
the span was only slightly increased. 

2. With slots in the sides, the variation of lift interference over the span was much 
less than that with the slots in the top and bottom walls. 

It should not necessarily be inferred from these results that for testing wings, slotted 
sides are preferable to slotted top and bottom walls. The variation of interference over 
the span must be balanced against the average value. If a choice of wind-tunnel geome- 
tries were available, it might be preferable to increase the width and decrease the height 
of the test section so that the tip of large-span wings would not be so near to the walls. 
To attain zero average lift interference in a square test section with solid top and bottom 
walls would require that the side walls be effectively almost open. 


As a check on the calculated interference for the large-span wing, an additional cal- 
culation was performed for a single lift element in the center of a square test section 
with slotted side walls. This lift element would correspond to a wing which spans only 
a very small portion of the tunnel width. The upwash interference factor for this case is 
given by 


a 





(29) 


This equation was derived independently but can be seen to be equivalent to equations (26) 
and (27) with x i? £j, y i? z i? and all equal to zero. Equation (29) gives for the 

upwash interference factor at the center 6 = 0.0469. This value is somewhat greater 
than that calculated near the center of the tunnel (near the midspan in fig. 4) for the large- 
span wing. Such an effect is to be expected because the interference factor for the large- 
span wing is a weighted average of contributions from lifting elements located along the 
span; and those nearer the tips are more affected by the slotted side walls than are those 
near the midspan position. Figure 5 shows how the upwash interference factor at the 
midspan position due to a lifting element falls off as the element is moved outward along 
the span. The interference at the midspan position due to an element located out along 
the span falls off still further when the wing is swept back. 
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of the single lift element 
h 


Figure 5-- Upwash interference factor at 
the center of the test section due to 
a single lift element in the horizontal 
center plane. 

The theory presented herein for the boundary-induced upwash interference due to 
lift on swept wings mounted in wind-tunnel test sections with slotted side walls and solid 
top and bottom walls is exact for the boundary conditions assumed. The accuracy of cal- 
culated results depends on accurate representation of the wing and on carrying the inte- 
grations and summations far enough to insure covergence. In consideration of the usual 
uncertainty in the knowledge of the lift distribution and of the approximate knowledge of 
the boundary conditions, it is believed that representation of the wing by ten lifting ele- 
ments as in the sample calculation is adequate. The use of homogeneous rather than 
exact slotted- wall boundary conditions is believed to introduce negligible error if there 
are several slots rather than only one or two in each slotted wall and if the wing is not too 
close to a wall. On the other hand, the action of the slots is uncertain and correspond- 
ingly .so is the true effective restriction constant. The restriction constant used is cal- 
culated on the assumption of potential flow outward from the test section through slots 
with thin edges and no boundary layer. For outflow through coarse slots, this method of 
calculation should yield approximately correct values for the restriction constant, but 
for very narrow slots or for strong inflow from the plenum chamber surrounding the slots, 
the calculated values of the restriction constant may be appreciably in error. If the slots 
are not of uniform width, an added uncertainty exists. 

In this paper, the computing equations are expressed in more general form than 
are those of reference 1; therefore, the model may be situated anywhere in the test sec- 
tion (provided it is not close enough to a boundary to invalidate the linearized theory 
assumed) and the upwash interference may be computed anywhere within the test section. 
The model must be so oriented that the lift vectors of the lifting elements lie approxi- 
mately parallel to the slotted side walls. 


t 
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Within the applicability of linearized theory, the upwash interference factor is not 
affected by compressibility at subsonic speeds, provided the stream boundary, including 
the slots, is approximately parallel to the tunnel stream direction. The theory presented 
herein may therefore be applied for subsonic compressible flow as well as for incom- 
pressible flow. 

The theory can be used to calculate the upwash velocities anywhere in the test sec- 
tion and thus corrections for the moment due to boundary-induced upwash at the tail or 
for the lift due to boundary-induced curvature of the flow can be applied. A compress- 
ibility correction can be made by applying the upwash velocities at x\l - rather than 
at x, where x is the distance from the lifting element and M is Mach number. The 
flow curvature must be multiplied by a compressibility factor — = For three - 

VI - M 2 

dimensional models of a practical size, the lift correction due to flow curvature is com- 
monly assumed to be negligible. 


RESUME 

The upwash interference factor due to wind-tunnel boundary interference on the lift 
of sweptback wings mounted in test sections with closed top and bottom walls and slotted 
side walls has been obtained in the form of infinite convergent integrals and summations 
suitable for calculation by means of a high-speed digital computer. The interference 
factor can be computed anywhere within the test section for an arbitrary wing mounted 
anywhere in the test section. In an example, the interference with slotted side walls is 
compared with that with slotted top and bottom walls. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 30, 1970. 
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^contribution to existing knowledge. . 


TECHNICAL TRANSLATIONS: Information 
published in aforeign language Considered 
to merit NASA distribution bn English, ✓ v 

SPECIAL PUBLICATIONS: Information 
derived from or of value ,to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks', 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION , , 
PUBLICATIONS: Information on technology 
used by ‘NASA that may be of particular : 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs* ^ 
Technology Utilization Reports and Notes, 
and Technology Surveys. 7 
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Details on the availability of these publications may be obtained from: 
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' SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 \ 


